Suppression of viscous fluid fingering: a piecewise constant-injection process 
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The injection of a fluid into another of larger viscosity in a Hele-Shaw cell usually results in the 
formation of highly branched patterns. Despite the richness of these structures, in many practical 
situations such convoluted shapes are quite undesirable. In this letter we propose an efficient and 
easily reproducible way to restrain these instabilities based on a simple piecewise constant pumping 
protocol. It results in a reduction in the size of the viscous fingers by one order of magnitude. 

PACS numbers: 47.15. gp, 47.20.Hw, 47.54.-r, 47.55. N- 



Hydrodynamic fingering is one of a wider class of in- 
terfacial instabilities occurring when one material is in- 
jected into another, or grows from a given chemical, bi- 
ological or geophysical process. Examples of such pat- 
tern forming systems include the classic Saff man- Taylor 
(ST) instability pL], the dynamics of chemically reacting 
fronts the growth of filamentary organisms [3], and 
lava flows These seemingly unrelated processes pro- 
duce complex patterns presenting similar morphological 
features, where highly branched shapes and dendriticlike 
configurations may arise. In this context, the control of 
the growth and form of the emerging complex morpholo- 
gies has long been a challenging topic, of great academic 
and technological relevance. 

Due to its relative simplicity, and multiple applications 
the ST instability has received much attention, and serves 
as a paradigm for pattern formation systems [5| . This in- 
stability arises when a fluid is injected against a second of 
much larger viscosity in the narrow gap between closely- 
spaced parallel glass plates of a Hele-Shaw (HS) cell [6|. 
Under constant flow injection rate, the result is the de- 
velopment of vastly ramified interfacial patterns (?- 
One of the most important practical situations related to 
this hydrodynamic instability is oil recovery |12|, where 
petroleum is displaced by injection of water into the oil 
field in an attempt to extract more oil from the well. In- 
terestingly, the dynamic behavior for flow in HS cells is 
described by the very same set of equations as those for 
flow in porous media [5]. In fact, the ST instability is a 
major source of poor oil recovery, once rapidly evolving 
fingers may reach the entrance of the well, and mainly 
water, and not oil is retrieved. This emblematic example 
clearly illustrates the importance of developing a funda- 
mental understanding of the dynamics of this type of sys- 
tem, and to find ways to contain, and possibly suppress 
such interfacial disturbances. 

Very recently, some research groups [iMi have ex- 
amined the possibility of avoiding the emergence of the 
usual branched morphology, by properly controlling the 
flow injection rate. Instead of considering a constant 
injection flux, they assumed a particular time- varying 
pumping rate which scaled with time like C(n)t~^/^, 
where C{n) depends on the interfacial wave number n. 



Their theoretical and experimental findings demonstrate 
that by using this specific pumping rate the formation 
of branched patterns is inhibited, and replaced by n-fold 
symmetric shapes. This process conveniently determines 
the number of emerging fingers. In spite of this con- 
trolling strategy, the resulting interfacial morphologies 
are notably noncircular, and still contain sizable fingers. 
Therefore, the interfacial fingering perturbations are not 
exactly wiped out, but redesigned into self-similar shapes 
with a prescribed number of fingers. In this sense, an 
efficient protocol for suppressing the development of the 
viscous fingering instability in radial HS flows is still lack- 
ing. 



Differently to what is done in Refs. [13l-ll5j. we con- 
sider a simple piecewise constant injection process which 
results in the actual suppression of the viscous fingering 
instability. In the usual constant injection situation, a 
certain amount of fluid is pumped in a finite time, and 
interfacial fingering results. In the procedure we suggest 
the average injection rate is kept unchanged, so that the 
same amount of fluid is injected in the same time interval, 
but interfacial irregularities are restrained. 

We begin by briefly describing the traditional radial 
flow setup in confined geometry. Consider a HS cell 
of gap spacing h containing two immiscible, incompress- 
ible, viscous fluids. The viscosities of the fluids are de- 
noted as rji and 7^2, and between them there exists a 
surface tension a. Fluid 1 is injected into fluid 2 at a 
constant injection rate Qo, equal to the area covered 
per unit time. Linear stability analysis of the prob- 
lem [^-[sl considers harmonic distortions of a nearly cir- 
cular fluid-fluid interface whose radius evolves according 
to IZiO^t) = R{t) + Cn(^) cos n6>, where the time depen- 
dent unperturbed radius is R{t) = Rt = ^ARo + Qo^/tt, 
represents the azimuthal angle, and n=0,±l, ±2, ... are 
discrete wave numbers. The unperturbed radius of the 
interface at t = is represented by i?o, and the Fourier 
perturbation amplitudes are given by 



Cn(t)=Cn(0)exp{/o(n)}, loin) = X{n)dt' (1) 
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where the Hnear growth rate is 



A(n) = 



(2) 



with /(n) = {A\n\ - l)/27r, g{n) = [b^a\n\{n^ - 
l)]/[12(7^i + 7^2)], and A = {r]2 - Vi)/{V2 + ^1) being the 
viscosity contrast. If Io{n) > the disturbance grows, 
indicating instabihty. In Eq. (j2j) we notice opposing ef- 
fects of the viscosity difference between the fluids (desta- 
bihzing) and of the surface tension (stabihzing). A rel- 
evant information can be extracted at the linear stage: 
the existence of a series of critical radii Rc{n) [or critical 
times tc{n)] at which the interface becomes unstable for a 
given mode n [defined by setting A(n) = 0], characteriz- 
ing a cascade of modes [9]. Therefore, to write the linear 
solution in a more realistic way, one should consider the 
interval [tc{'n)^tf] only, because integration over [0,tc(n)] 
would lead to an artificial diminishing in the size of the 
fingers, so that Cn(^c(^)) < Cn(0). This is an unphysical 
effect because we assume that the amplitudes can not go 
below Cn(0) due to noncontrollable factors, e.g., irregu- 
larities on the surface of the plates ^] . For this reason we 
will consider the initial perturbation to be independent 
of n, that is, Cn(0) = Co- Notice, however, that for lower 
modes tc{n) <^tf^ and solution ([T]) can be used as a good 
approximation. 

We proceed by describing the stabilization protocol. 
Consider a radial HS flow, and suppose that one is re- 
quired to pump a fluid into the cell, at a specified av- 
erage rate Qo, during the time interval [0,t/]. This 
defines the final area occupied by the injected fluid, 
Af = ttRq + Qotf. Given the quantities Af and t/, 
our goal is to design a simple injection process that sup- 
presses viscous fingering events occurring at usual con- 
stant pumping procedure which utilizes the same input 
parameters Af and tf. The results presented here refer 
to Af ^ 140 cm^ and t/ = 28.0 s (Qo = 5.00 cmVs) and 
to the following set of characteristic physical parameters: 
Ro 0.30 cm, b 0.10 cm, Co i^o/2400 cx lO""^ cm, 
a = 63.0 dyne/cm, r]2 = 5.21 g/cm s, and r]i ^ 0. These 
values are consistent with those used in typical experi- 
mental realizations We chose A ^ 1 because in 
this limit the most unstable situation is reached, so that 
the results can only improve for < A < 1. 

Out of a plethora of possibilities for time- varying injec- 
tion scenarios, the simplest nontrivial procedure is a two- 
stage piecewise constant pumping (Fig. [1]), whose choice 
will be justified shortly. Specifically, we split the time 
interval [0,t/] into [0,r], during which the constant in- 
jection rate is Q; and [r, t/], with tf = (l + /3)r, for which 
the injection is given by 7(5. The requirement that the 
average injection remains unchanged demands that 
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FIG. 1. (Color online) Injection rate Q(t) as a function of time 
for a two-stage piecewise constant protocol. The equivalent 
constant injection rate Qo is represented by the horizontal 
dashed line. 



with P > and 7/3 > —1. If the parameter 7 is negative 
we have an injection stage followed by a period of suction. 
For 7 > two nontrivial possibilities arise: if 7 < 1 the 
injection in the first stage is stronger than in the second, 
and for 7 > 1 the weaker injection stage comes first. For 
7 = 1 we recover the usual constant pumping situation. 

As commented earlier, the formal output of a linear 
analysis is given by Eq. ^ for the constant injection, 
while for the lower modes in the two-stage pumping we 
have Cn(^) — Coexp{/'(n)} with 

r{n) = Ai(n)dt + J \2{n)dt = 27Tf{n)A{Rf,Ro) 



2^^(n)(l + /37) 
Qo(l + /3) 



-T{Rf,Rr)^T{Rr,Ro) 



(4) 



(1 + 7/5) 



Qo 



(3) 



where Ai and A2 refer to the first and second stages with 
their respective injection rates, A{x^y) = ln{x/y), and 
r{x,y) = {1/y — 1/x). Note that while Rq and i?/, 
the initial and final unperturbed radii, are constant, the 
radius at t = r, is a function of the free parameters, 

Rr = Rr{i3,j). 

If we are to suppress instabilities, at least for the initial 
modes, we must impose that Cn/Cn = exp{/' — /q} < 1 in 
the end of the whole process, i.e., at t = t/. By analyzing 
the overall sign in the argument of the exponential one 
shows that the only scenario that leads to stabilization 
is that with 7 > 1, or, an initial stage with a relatively 
weak injection rate followed by a stronger (by a factor of 
7) injection stage. This can be understood by recalling 
that the lower the mode the larger its maximum growth 
rate. So, to suppress the initial and larger instabilities 
a slow injection is needed. For the other scenarios no 
values of f3 and 7 produce a stabilizing effect. We have 
also checked oscillating pumping with Q{t) > and os- 
cillations involving injection and suction in each cycle, 
both resulting in an enhancement of the ST instability. 
Thus, our task is to establish a general and simple way 
to obtain optimal values, /3* and 7*, within the selected 
scenario (/3 > and 7 > 1). 

It is a key point to realize that during the piece- 
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wise constant injection the cascade of modes occurs dif- 
ferently. In general, the "wakening" of lower modes 
still occurs in the first stage of pumping. However, for 
t = tf/{l -\- P) = r the injection rate jumps discontin- 
uously and two effects occur: the most unstable modes 
prior to the sudden change in the pumping regime have 
their growth rates drastically decreased, and, at the same 
time, a certain number of new modes become unstable in- 
stantaneously. After that the cascade process continues. 
We have used a smoothed out version of the step-like in- 
jection function depicted in Fig. [T] and verified that the 
outcome is virtually unchanged. This signifies that our 
results are not due to any pecuHarity related to the dis- 
continuous nature of Q{t). 

In our quest for the optimal parameters /3* and 7*, we 
try to devise simple rules that are independent of the dis- 
crete wave number. For this reason, we will suppose that 
the relevant modes approximately fit in one out of two 
classes: (i) modes that become unstable soon after the 
beginning of the first stage (lower modes), and (ii) modes 
which attain their regimes of instability exactly at t = r 
or a bit later (higher modes). The remaining modes, i.e., 
those which become unstable at the end of first stage 
[t ^ tf /{I -\- P)] and in the final part of the second stage 
{t ^ tf) are not critical because their regime of strong in- 
stability lasts for a short time. The effectiveness of these 
considerations will be made evident in what follows. For 
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FIG. 2. (Color online) Behavior of I'^a: 
of 7. The absolute minimum at 7* = 3.72 is indicated. 

the relevant modes of type (i) we have the perturbations 
appropriately described by Eq. (j4|), while for the modes 
of type (ii) we have Cm(t) = Co exp{/''(m)}, where 



J"(m) 



\2{m)dt = 27r/(m)A(i?/,i?o) 
27r<?(m)(l + /?7), 



7Qo(l 



-r{Rf,Rr) 



(5) 



In order to get rid of the dependence on n and m, we will 
focus on the modes associated to the largest perturbation 
amplitudes in each stage, Umax and rumaxi respectively. 
Therefore, we can guarantee that the size of the fingers 
corresponding to these wave numbers constitute an upper 



bound for the scale of all other modes. This is obtained 
by setting dV jdn — and dl" /dm = 0, yielding 



'1 



2Qo(??i+r?2)7(l + /3)A(i?/,i?o) 



3 7r62CT(l + /37) [r(i?/, Rr) + iT{Rr, i?o)] ' 

(6) 



and 



/I ^ 2Qo(r/i + 7/2)7(1 + /9)A(%, Rr) 



(7) 



7r62a(l + /37)r(i?/,i?,) 
The referred upper bounds are then I'{nmax) and 
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FIG. 3. (Color online) The continuous curves represent I' and 
I" as functions of n, calculated for /3* and 7*. The dashed 
curve shows the corresponding quantity /o for the constant 
injection process. 

I"{ramax)- Before imposing any minimization condition 
on these quantities we make sure that the suppression 
of the instabilities is uniform, with higher modes as con- 
trolled as the lower ones, by setting the constraint 

which reduces the dimension of our space of parameters 
and can be numerically solved to give /3(7). Using this 
relation to eliminate f3 in /^^^(/3,7), or I!^^^{f3,j), we 
obtain a function of the single variable 7, shown in Fig.[2l 
Finally, we pick the value of this variable that minimizes 
^max [^^^ automatically I'^^^x ®]' ^^^^ corre- 

sponds to the optimal value 7* = 3.72 and also yields 
P* = P{j*) = 0.46. This completely characterizes the 
piecewise constant injection and gives r ^ 19.2 s be- 
fore which the injection rate is Q = 2.7 cm^/s, with the 
stronger pumping lasting for about 8.8 s with an injec- 
tion of 7(5 = 10 cm^/s. In Fig. [3] we show /'(n), /''(n) 
(optimal parameters used) , and /q for the equivalent con- 
stant pumping process as functions of the wave number. 
Condition (|8]) is evident from this figure and, most im- 
portantly, we have lomax ^ 7.5 and I'^^^ 1'^^^ ^ 4.9. 
Since these integrals are related to the logarithm of the 
amplitude, the decrease in the relative size of the largest 
fingers is of one order of magnitude. The effectiveness 
of the protocol can be seen even more clearly in Fig. |H 
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FIG. 4. (Color online) Comparison between the interfacial 
patterns formed during constant injection (a), and piecewise 
constant injection (b) including 40 Fourier modes and a ran- 
dom choice of phases. The same is valid for figures (c) and 
(d), respectively, but considering 30 modes, and a different 
set of random phases. 



where we plot the linear evolution of the interfaces for 
the equivalent constant pumping in (a) and (c), and for 
the two-stage pumping in (b) and (d). The patterns on 
the top have the same initial conditions (including the 
random phases attributed to each mode [8, 9]), and 40 
Fourier modes have been considered. The same is valid 
for the two bottom panels, but including 30 modes and 
a distinct set of random phases. The interfaces are plot- 
ted in intervals of t//5. Note that they are spaced in a 
more uniform way for the constant injection, while for the 
piecewise constant injection the interfaces are initially 
closely spaced, and then more widely separated. We em- 
phasize that the cascade of modes was considered with- 
out any simplifying hypothesis in these numerical results. 
Moreover, notice that the shapes shown in Fig. IH^b) and 
Fig. m^d) are very similar, revealing an insensitivity to 
changes in the initial conditions. 

We have also investigated the weakly nonlinear evo- 
lution of the interfaces by considering both second and 
third order couplings [9], and verified that our protocol 
produces patterns nearly identical to those depicted in 
Fig. m^b) and Fig. ID^d). This indicates that the emer- 



gence of nonlinearities is unfavored, expanding the dura- 
tion of the linear regime. However, the robustness of this 
stabilization for fully nonlinear stages of the dynamics 
merits further numerical and experimental investigation. 

In conclusion, we have introduced a simple injection 
process for which interfacial viscous fingering instabili- 
ties are truly suppressed. The procedure does not rely 
on unusual material properties of fluids, or on drastic 
modifications of the traditional radial HS flow setup. It 
only requires the employment of an optimal two stage 
piecewise constant injection mechanism. This stabiliza- 
tion strategy might be useful to improve the efficiency 
and control of a number of physical, biological, and tech- 
nological problems related to viscous fingering phenom- 
ena. 
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